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Abstract. We theoretically study the collective excitations of an ideal gas confined in an isotropic harmonic 
trap. We give an exact solution to the Boltzmann-Vlasov equation; as expected for a single-component 
system, the associated mode frequencies are integer multiples of the trapping frequency. We show that the 
expressions found by the scaling ansatz method are a special case of our solution. Our findings, however, 
are most useful in case the trap contains more than one phase: we demonstrate how to obtain the oscillation 
frequencies in case an interface is present between the ideal gas and a different phase. 
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1 Introduction 

Since the experimental realization of quantum degenerate 
gases, the role of particle interactions has been thoroughly 
studied. While weak interactions characterize most ultra- 
cold systems due to their diluteness, the use of Feshbach 
resonances allows the study of very strongly interacting 
systems. Although interactions are generally important, 
many experimentally realized gases may be well-modelled 
by ideal gases. For example, collisions are suppressed in 
fully-polarized Fermi gases at temperatures well below 
their Fermi temperature because s-wave interactions are 
forbidden due to Pauli blocking. Such gases appear for in- 
stance in the experiments on imbalanced fermion gases [H 
El[3] and in phase-segregated Bose-Fermi mixtures [2J[5l[H] • 
A description in terms of ideal gases is also appropriate 
in ultracold two-component Fermi systems when the in- 
terspecies scattering length is tuned to zero by means of 
a Feshbach resonance [7]. Moreover, it is well-known that 
sufficiently dilute gases become collisionless; as a result, 
collisionless behavior also sets in as the temperature of 
trapped gases well above the temperature of quantum de- 
generacy [8]. 

For a hydrodynamic gas (or superfluid) under exper- 
imentally relevant conditions, the variables relevant for 
studying the collective excitations are, for example, the 
local density and velocity; on the other hand, for a col- 
lisionless fluid, one must, in general, find the entire dis- 
tribution function in momentum- and coordinate-space. 
In some cases, one may extract dynamical solutions us- 
ing a scaled form of the equilibrium distribution func- 
tion [9lfl ^ fTT 1 [T2 l lT3] . However, the scaling approach fails 
for a two-component system with an interface • Other 
methods to obtain the collective mode frequencies include 



the method of averaging, the sum-rule approach, and the 
random-phase approximation p^[T^[T6lJT71fT8l[T9| [2C7| . 

We present here a solution for the Boltzmann-Vlasov 
equation of an ideal gas in an isotropic harmonic trap. 
For the single-component system, we recover the spec- 
trum nujQ with n integer and loq the trapping frequency. 
We then show how our solution may be used to study a 
two-component system with an interface by writing down 
the appropriate boundary conditions at the interface. For 
the monopole modes, we give analytical expressions which 
may be straightforwardly used to calculate their frequen- 
cies in isotropic and highly-elongated elliptical traps. Fi- 
nally, we discuss damping of these modes. 



2 Boltzmann-Vlasov Equation 

An ideal gas, confined by an isotropic potential tow 2 ! -2 /2, 
is described by a distribution function /(r,p,i), evolving 
according to the Boltzmann-Vlasov equation: 



md t f + p • V P / - m^r • V p / = 



(1) 



In the following fo denotes the equilibrium solution, for 
which dtfo = 0. Collective oscillations may be studied 
by looking at periodic deviations from /q. Writing / = 



fo 



we obtain the result 



■■C{{L h E s y) 



a 3 (j> 3 



o/2luq 



(2) 



with the a.j arbitrary dimcnsionless constants and C an 
arbitrary function of the angular momentum components 
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Lj = (r x p)j and of the Ej = (p 2 /m + muj 2 r 2 )/2 (with 
j = x,y,z). Note that solution ([2]) is not obtained by 
linearization so v may be large compared to fo- 



3 Single-component System 

Consider a trap containing a single particle component. 
Due to the absence of boundary conditions, the function 
C in Eq. ([2]) is arbitrary; physically, v/C must then be con- 
tinuous everywhere. Hence, lo has to be chosen such that 
no discontinuities appear when the coordinates Tj cross 
zero. This is guaranteed if lo = uloq (from now on, n is 
an integer) , which is the well-known spectrum of a single- 
component trapped collisionlcss gas 8j . We show now that 
indeed all experimentally relevant collective modes can be 
obtained from Eq. @ . 

The lowest dipole or Kohn mode can be found by tak- 
ing, for example, a = (0, 0, 1); it is then clear that lo = loo 
ensures that vjC converges to the same value if r z — > + 
or r z — > 0~ when p z < 0. It is also easy to see that, when 
C is rotationally invariant in coordinate space, the den- 
sity deviation Sp oc J d 3 p^ is simply proportional to r z , 
which, as expected for the dipole mode, is proportional to 
the spherical harmonic Y^_ x . 

Due to rotational symmetry in space the monopole 
modes can be described by taking a = (1, 1, 1) such that 

v = C (L 2 , E) [p 2 - m 2 c; 2 r 2 + 2m W „r • p] 1 "'*" , 

with E = J2j Ej ■ Arguments similar to those for the 
dipole mode lead to the conclusion that lo — 2uloq, which 
is the well-known spectrum for monopole modes. More- 
over, for the lowest excitation (n — 1), the density devia- 
tion Sp is proportional to the spherical harmonic Y®_ . 

Finally, for the quadrupole mode, two of the dj must 
be equal while the third opposite in sign. Choosing a = 
(1,1, —2), one finds that lo = 2w and that Sp oc x 2 + y 2 — 
2z 2 oc Y£_ 2 , as it mustQ. 

A highly-elongated elliptical trap may be modelled 
by an isotropic 2D confinement U(r) — mcu 2 (r 2 + r 2 )/2. 
This results in a cylindrical configuration with the a;-axis 
along the axial direction and rotational symmetry in the 
y — z plane. The angular dependencies of Sp for the low- 
est monopole, dipole and quadrupole modes are known to 
vary as e M oc (y + iz) e with 9 the angle in the y — z plane 
and £ = 0, 1, 2. The Kohn mode (£ = 1) can therefore be 
recovered by an appropriate superposition of Eq. with 
a = (0, 1, 0) and Eq. © with a = (0, 0, 1), together with a 
rotational-invariant function C. In an analogous manner, 
and using a C-function which depends on L Xl one may 
obtain the quadrupole mode {1 = 2). 

We now explain the relation between Eq. (|2|) and the 
scaling ansatz method, widely used in the literature [9. 10. 

EHE]. 

1 Note that for this to be true, it is necessary that C is rota- 
tionally symmetric. 



4 Relation to the Scaling Ansatz 

Consider a confined single-component gas with equilib- 
rium distribution /q. Due to the radial symmetry, /q can 
be written as a function of L 2 and E only; we assume that 
fo is a function of E only. According to the scaling ansatz, 
the distribution function / of a trapped single-component 
system at a certain time t may be written in terms of the 
equilibrium distribution fo as follows: 

/(r,p,*)=/o(r(t),p(t)). (3) 

In order to find the monopole and the quadrupole modes 
one continues by takes the scaling functions rj(t) = rj/ctj 
and Pj(t) = cijPj — mctjrj where the dot denotes the 
derivative with respect to time and otj is a periodic func- 
tion of time [9 , 1 Ol iTTUTS] . For small deviations from equi- 
librium, one writes aj = 1 + eaje~ luJt with e <C 1 and aj 
dimcnsionless constants. Linearization of / in terms of e 
yields: 

f = fo+ee-^ t m- 1 (d E f ) (4) 

x I [ a j(Pj + iniu>o r j) 2 + imaj(uj — 2uio)rjPj\ 

\j=x,y,z j 

Clearly / reduces to Eq. ([2]) when lo = 2lvo, which, as 
we have seen, is indeed the frequency associated with the 
lowest monopole and quadrupole modes of an ideal gas. 

The Kohn mode, on the other hand, may be recovered 
by substitution of r z (t) — r z — (3 and p z (t) = p z — P while 
for j = x, y, one takes rj(t) — rj and pj(t) — pj. Assuming 
a time-dependent function (3{t) — ee~ luJt and linearizing 
in terms of e yields: 

/ = f + iee- wt {d E fo) [up z + imij 2 r z ] . (5) 

Again this reduces to our solution ([2]) in case a = (0, 0, 1), 
and lo = ci>q. 

We therefore conclude that, in a linearized form, the 
expressions for the distribution functions used by the scal- 
ing ansatz method reduce to our solution ([2]) in case lo = 
loo and lo = 2loq. Note that, while our solution (J5J) solves 
the Boltzmann-Vlasov equation, Eqs. (j4|) and |5| do not, 
except for specific values of lo. Therefore, the reason why, 
for a single-component system, the correct spectrum arises 
from Eq. (|2|), and from Eqs. (J4j) and (O, is different. Note 
further that, although the scaling ansatz does not give the 
most general solution for the Boltzmann-Vlasov equation, 
it may be generalized to find oscillation frequencies using 
a Boltzmann equation which includes a mean-field and 
collision term [9 | fT0 U TT V T2] . 

5 Two-Component System with an Interface 

Assume in the following the presence of two phases, sepa- 
rated in the trap by an interface at radial position £. We 
investigate the boundary conditions at the interface for 
the ideal-gas phase and leave the nature of the other phase 
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unspecified. Experimentally relevant cases include the im- 
balanced fermion gases [Hf2lf3l[2T| and phase-segregated 
Bose- Fermi mixtures [4,5,6,22,23]. 

Under the assumption that particles are specularly re- 
flected, conservation of energy and particle number relate 
the distribution function of incoming and outgoing parti- 
cles at the interface |24) : 

KC,x)-KC J -x) = 2mxC(^/o), (6) 

with x defined by r • p = x r P where r = |r| and p = |p|. 

Consider the point r = (0, 0, () on the interface. The 
sole result of a particle reflecting there is its change of ve- 
locity component p z to — p z . Therefore, aside from L z , E x 
and E y , also L x , L y and E z are invariant under reflection, 
such that C(x) = C(~x)- The solution of Eq. © satisfy- 
ing the boundary condition (|6|) at position r = (0, 0, £) is 
then: 

V = 2mxt(d E h)[l-e-^]-\ (7) 

where, 

t = uJ 1 Arg [a x pl + a y p 2 y + a z (xp + zmw C) 2 ] ■ (8) 

A general solution for v at r = (0, 0, Q can be obtained 
by taking superpositions of Eq. ([7]) with each different 
values of a. In the following, we show how one can study 
breathing modes for gases confined in isotropic and highly- 
elongated harmonic traps. Multipole modes may then be 
studied in an analogous manner. 

Spherical trap — The breathing modes can be found by 
imposing the condition of rotational symmetry in Eq. ([8]) , 
such that: 

r = uJq 1 Arg [p 2 - m 2 uj 2 C 2 + 2imLj a x(p] ■ (9) 

Physically, r is the time for a particle departing radially 
from the interface, to return to it. For a fully-degenerate 
Fermi gas, the values for £ and lo for which damping occurs 
are depicted in Fig. 2 of Ref. [6]. Note that, due to the 
rotational invariance Eq. ([7|) is valid at any point on the 
interface. 

Highly Elongated Traps — In case of highly-elongated 
harmonic traps, one can approximate the trapping poten- 
tial by U(r) — mujo(r 2 + r 2 )/2. The breathing modes can 
be studied assuming a = (1, 1, 0) and considering v in the 
spatial point (0, 0, £) where Eq. ([7]) is valid with: 

r = Wo x Arg [p 2 ( X 2 + (1 - X 2 ) cos 2 <j>) 

-m 2 ulC 2 + 2imLu oX Cp] ■ (10) 

Here <p is the azimuth of p in the x-y-z coordinate system. 

Apart from the boundary condition ([6]), which by it- 
self ensures that the mean velocity of the coUisionless gas 
at the interface is equal to another boundary condi- 
tion applies, ensuring local mechanical equilibrium at the 
interface. This condition is the Laplace equation and it 
relates the radial pressure tensors U rr of the inner phase 
A and the outer phase B at the interface: 

^-77,^=7^(1/^+1/^). (11) 



Here R t and R 2 are the principal radii of curvature of the 
interface and jab is the interface tension of phase A and 
B. 

Knowledge of the distribution function v from Eq. (J5J) 
allows a direct calculation of the pressure tensor TI rr in 
Eq. (fTT|) for the ideal-gas phase. Assuming one knows the 
pressure tensor for the other phase, all the necessary ingre- 
dients are present to find the collective mode frequencies, 
as shown in Refs. [3ll21j. 

6 Work Done by a Moving Interface 

We now show that the total energy transfer through the 
interface over one oscillation period vanishes unless a pole 
is present in the pressure tensor, in which case the collec- 
tive mode gets damped. We determine also a criterion for 
this damping to occur. 

The work done by a moving interface onto a collision- 
less gas is most simply expressed using the pressure tensor 
normal to that interface, EI rri and the average velocity of 
the gas in the same direction, u r (here r denotes the ra- 
dial coordinate). At r = (0,0, C), the functions Sn rr and 
u r vary as: 

6II rr ex e-^J^dpj d X xV(a B /o)cot (^) 

and 

poo 

u r ^ie- Mt dpp 3 (d E f ), 
Jo 

where we have used that t(x) = ~ T (~ x)- The work done 
per unit area during one oscillation period T = 27t/K(oj) 
is given by 

aw = [ dm(sn rr )u(u r ). 

Jo 

If there appears no pole in the integrand of 6n rr , then 
dn rr and u r are 7r out of phase and the total work done 
during a single cycle averages to zero. If, on the other 
hand, the integrand of 5n rr has a pole for some x, then 
one must pass below that pole when integrating over x [25] ; 
as a result, 5n rr acquires an imaginary part and is no 
longer n out of phase with u r . This results in net work 
being done on the coUisionless gas during one period. In 
other words, it corresponds to damping of the collective 
modes, with their energy transferred to single-particle modes 
in the coUisionless gas (in analogy to Landau damping) . 

One may ask now when a pole appears in the integrand 
of 5n rr . From Eq. ([7]), ^(C,x) diverges if: 

r(C,x)=27rn/w, (12) 

with n again an integer. This gives rise to the following 
criterion: damping at the interface can only happen when 
the collective mode frequency exceeds 2luo. Indeed, since 
r| < tt/ujq, Eq. (| 1 2[) can only be satisfied when lo > 2wo 

Finally, note that it can be shown that, if no pole ap- 
pears in the integrand of Sn rr at the interface, then there 
is no pole anywhere in the trap. 
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7 Conclusions 

We have presented a solution (see Eq. ©) for the Boltz- 
mann-Vlasov equation, describing the collective excita- 
tions of an ideal gas in an isotropic harmonic trap. If the 
trap is entirely filled by the ideal gas, the associated mode 
frequencies are simply integer multiples of the trapping 
frequency. In that case, we proved that the expressions 
for the distribution function used by the scaling ansatz 
method, are a special case of our solution. We further 
applied our solution to a trap consisting of an ideal gas 
in contact with a different phase by means of an inter- 
face. Taking the particles of the ideal-gas phase to specu- 
larly reflect on the interface, we show how the breathing 
mode frequencies may be obtained in isotropic and highly- 
elongated elliptical traps. Finally, we discuss how, in the 
presence of an interface, damping may arise. 
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